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LONG-TERM  GOALS 

We  are  part  of  a  multi-institutional  researeh  team*  funded  by  the  ONR-sponsored  Radianee  in  a 
Dynamie  Ocean  (RaDyO)  program.  The  primary  research  goals  of  the  program  are  to  (1)  examine 
time-dependent  oeeanic  radianee  distribution  in  relation  to  dynamie  surfaee  boundary  layer  (SBL) 
proeesses;  (2)  eonstruet  a  radianee-based  SBL  model;  (3)  validate  the  model  with  field  observations; 
and  (4)  investigate  the  feasibility  of  inverting  the  model  to  yield  SBL  conditions.  The  goals  of  our 
team  are  to  eontribute  innovative  measurements,  analyses  and  models  of  the  sea  surface  roughness  at 
length  seales  as  small  as  a  millimeter.  This  eharaeterization  ineludes  mieroseale  and  whiteeap  breaking 
waves. 
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OBJECTIVES 

Nonlinear  interfaeial  roughness  elements  -  sharp  erested  waves,  breaking  waves  as  well  as  the  foam, 
subsurface  bubbles  and  spray  they  produee,  eontribute  substantially  to  the  distortion  of  the  optical 
transmission  through  the  air-sea  interfaee.  These  eommon  surfaee  roughness  features  oeeur  on  a  wide 
range  of  length  seales,  from  the  dominant  sea  state  down  to  eapillary  waves.  Wave  breaking  signatures 
range  from  large  whiteeaps  with  their  residual  passive  foam,  down  to  the  ubiquitous  eentimeter  seale 
mieroseale  breakers  that  do  not  entrain  air.  There  is  substantial  complexity  in  the  loeal  wind-driven  sea 
surfaee  roughness  mierostrueture.  Traditional  deseriptors  of  sea  surfaee  roughness  are  seale-integrated 
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statistical  properties,  such  as  significant  wave  height,  mean  squared  slope  (eg.  Cox  and  Munk,  1954) 
and  breaking  probability  (e.g.  Holthuijsen  and  Herbers,  1986).  Subsequently,  speetral  eharaeterisations 
of  wave  height,  slope  and  eurvature  have  been  measured,  providing  a  seale  resolution  into  Fourier 
modes  for  these  geometrieal  sea  roughness  parameters.  More  reeently,  measurements  of  whiteeap  erest 
length  spectral  density  (eg.  Phillips  et  al,  2001,  Gemmrieh,  2005)  and  mieroseale  breaker  erest  length 
speetral  density  (eg.  Jessup  and  Phadnis,  2005)  have  been  reported. 

Our  effort  seeks  to  provide  a  more  eomprehensive  deseription  of  the  physieal  and  optieal  roughness  of 
the  sea  surface.  We  will  aehieve  this  by  implementing  a  eomprehensive  sea  surface  roughness 
observational  ‘module’  within  the  RADYO  field  program  to  provide  optimal  eoverage  of  the 
fundamental  optieal  distortion  proeesses  assoeiated  with  the  air-sea  interfaee.  Within  our  innovative 
eomplementary  data  gathering,  analysis  and  modeling  effort,  we  will  pursue  both  speetral  and  phase- 
resolved  perspectives.  These  will  eontribute  directly  towards  refining  the  representation  of  surface 
wave  distortion  in  present  air-sea  interfaeial  optieal  transmission  models. 

APPROACH 

We  will  build  substantially  on  our  aeeumulated  expertise  in  sea  surfaee  proeesses  and  air-sea 
interaetion.  We  are  working  within  the  larger  team  (listed  above)  measuring  and  eharaeterizing  the 
surface  roughness.  The  group  plans  to  eontribute  the  following  eomponents  to  the  primary  sea  surfaee 
roughness  data  gathering  effort  in  RaDyO: 

•  polarization  camera  measurements  of  the  sea  surfaee  slope  topography,  down  to  eapillary  wave 
seales,  of  an  approximately  Im  x  Im  pateh  of  the  sea  surfaee  (see  Figure  1),  captured  at  video  rates. 
[Sehultz] 

•  co-located  and  synchronous  orthogonal  75  Hz  linear  scanning  laser  altimeter  data  to  provide  spatio- 
temporal  properties  of  the  wave  height  field  (resolved  to  0(0. 5m)  wavelengths)  [Banner,  Morison] 

•  high  resolution  video  imagery  to  reeord  whiteeap  data,  from  two  eameras,  elose  range  and  broad 
field  [Gemmrieh] 

•  fast  response,  infrared  imagery  to  quantify  properties  of  the  mieroseale  breakers,  and  surfaee  layer 
kinematies  and  vortieity  [Zappa] 

•  sonic  anemometer  to  characterize  the  near-surface  wind  speed  and  wind  stress  [Zappa] 

Our  envisaged  data  analysis  effort  will  inelude;  detailed  analyses  of  the  slope  field  topography;  laser 
altimeter  wave  height  and  large  seale  wave  slope  data;  statistieal  distributions  of  whiteeap  erest  length 
density  in  different  seale  bands  of  propagation  speed  and  similarly  for  the  mieroseale  breakers,  as 
funetions  of  the  wind  speed/stress  and  the  underlying  dominant  sea  state.  Our  eontributions  to  the 
modeling  effort  will  foeus  on  using  the  data  to  refine  the  sea  surfaee  roughness  transfer  function.  This 
comprises  the  representation  of  nonlinearity  and  breaking  surfaee  wave  effects  including  bubbles, 
passive  foam,  aetive  whiteeap  eover  and  spray,  as  well  as  mieroseale  breakers. 
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WORK  COMPLETED 


Our  effort  in  FY07  has  been  primarily  in  the  detailed  planning  of  the  suite  of  sea  surface  roughness 
measurements  that  we  will  undertake  during  the  Scripps  Institution  of  Oceanography  (SIO)  Pier 
Experiment  scheduled  for  January  6-28,  2008.  During  FY07  we  refined  our  choices  of  the 
instrumentation  needed  to  make  the  measurements  described  in  the  preceding  section,  and  continued 
work  on  the  analysis  techniques  for  characterizing  the  various  roughness  features.  We  participated  in 
the  FY07  RaDyO  scientific  planning  meetings,  which  were  held  in  Montreal  in  October  2006  and  at 
SIO  in  June  2007. 

We  carried  out  data  processing  and  validation  of  our  scanning  lidar  data  gathered  at  the  USAGE  Duck 
pier  field  site  at  Duck,  N.  Carolina.  Two  of  these  instruments,  operating  in  quadrature,  will  be 
deployed  in  the  RaDyO  field  experiments  to  monitor  the  large  scale  wave  geometry  (height  and  slope 
components).  We  also  progressed  with  our  development  of  a  robust  ‘individual  wave’  decomposition 
capability  so  that  local  physical  roughness  elements  can  be  detected  and  characterized  along  with  their 
space-time  phasing,  thereby  overcoming  the  classical  Fourier  spectrum  issue  of  bound  versus  free 
wave  contributions  in  assessing  true  physical  sea  surface  roughness. 

RESULTS 


(i)  Instrumentation  proposed  for  the  SIO  pier  experiment,  January  6-28,  2008. 

Figure  I  below  shows  schematically  the  instrumentation  to  be  deployed  in  this  field  testing  phase. 

2008  SIO  PIER  EXPERIMENT  INSTRUMENTATION  SCHEMATIC 
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video  _ 
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PLAN  VIEW 


orthogonal  scanning  Ildars  Infrared  camera 


Figure  1.  Schematic  of  instrumentation  packages  to  be  deployed  from  the  northern  swinging  boom 
at  the  end  of  the  Scripps  pier.  The  end  of  the  instrumentation  boom  will  be  about  9m  from  the  edge 
of  the  pier  and  about  10m  above  the  mean  water  level  The  approximate  field  of  view  of  the  various 
instruments  is  shown.  Another  wide  angle  whitecap  video  camera  is  mounted  well  above  the  pier. 
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Banner/Morison  plan  to  deploy  two  orthogonal  line  scanning  Ildars,  synchronized  for  zero  crosstalk. 
These  will  be  positioned  on  the  boom  so  that  their  intersection  point  is  within  the  common  footprint  of 
the  polarimetric  (Schultz),  infrared  (Zappa)  and  visible  (Gemmrich)  imagery  cameras  to  measure 
small-scale  surface  roughness  features  and  breaking  waves. 

Zappa  will  deploy  his  infrared/visible  camera  system  (with  blackbody  target,  a  blackbody  controller,  a 
laser  altimeter).  He  will  also  deploy  his  environmental  monitoring  system  (sonic  anemometer,  a  Licor 
water  vapor  sensor,  a  Vaisala  RH/T/P  probe,  a  motion  package,  a  pyranometer,  and  a  pyrgeometer). 

Gemmrich  will  deploy  2  video  visible  imagery  cameras.  One  camera  will  be  mounted  on  the  main 
boom  next  to  our  other  instrumentation  packages.  The  second  camera  will  be  mounted  higher  up  to 
provide  a  wider  perspective  on  larger  scale  breaking  events. 

Schultz/Corrada-Emmanuel  will  deploy  an  instrument  package  located  on  the  boom  that  includes  a 
polarimetric  camera  imaging  the  very  small-scale  waves,  an  autofocus  controller  for  this  camera,  a 
laser  rangefinder  for  the  autofocus  mechanism,  a  polarimetric  camera  looking  up  at  the  sky  and  a 
motion  package. 

The  individual  data  acquisition  systems  will  be  synchronized  to  GPS  accuracy  so  that  the  various  data 
sets  can  be  interrelated. 

(ii)  Duck  pier  scanning  lidar  breaking  wave  profiles 

Our  scanning  lidar  system  was  field-tested  at  the  Duck  pier  site  over  a  wide  range  of  conditions  where 
the  wind  speed  Uio  ranged  from  light  and  variable,  up  to  16  m/s.  Figure  2  shows  typical  scanning  lidar 
data  over  breaking  waves  under  strong  wind  forcing  conditions. 
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(a) 


(c) 


(e) 


(f) 


Figure  2.  Panels  (a)  through  (f)  show  the  temporal  evolution  of  the  surface  elevation  profile 
measured  by  our  scanning  lidar  of  a  breaking  wave  shoaling  in  7m  water  depth  at  the  Duck  pier 
site.  The  observed  wave  geometry  shows  a  very  steep  leading  face,  with  a  maximum  steepness 
approaching  45°.  Note  the  persistent  feature  developing  on  the  forward  face.  The  x  and  z  axes  show 
distances  in  meters.  The  time  interval  between  the  successive  images  is  0.1 76  seconds. 


The  lidar  was  deployed  from  a  boom  at  a  height  of  9m  above  the  mean  sea  level,  where  the  water  depth 
was  about  7m.  The  lidar  sensing  spot  on  the  sea  surface  had  a  diameter  of  ~0.2m,  it  made  an  altimetric 
height  determination  every  0.5°  and  the  temporal  scan  rate  was  in  excess  of  50  Hz.  With  its  viewing 
window  of  ±45°  to  the  vertical,  the  maximum  horizontal  data  extent  was  ±  9m  from  nadir. 

We  deployed  this  instrument  for  several  days.  As  anticipated,  the  horizontal  extent  of  the  lidar  return 
increased  as  the  wind  strengthened,  because  of  the  greater  number  of  scatterers  that  provide  the 
specular  return  facets  for  the  lidar  to  make  its  time-of-flight  measurement.  Our  experience  confirmed 
that  this  method  will  provide  useful  data  on  the  height  and  local  directional  slope  of  the  dominant 
waves.  This  information  characterizes  the  background  environment  experienced  by  the  short  wind 
waves  (the  sea  surface  microstructure  roughness).  This  information  will  also  allow  accurate  phasing  of 
the  polarimetric  camera  imagery  of  the  sea  surface  microstructure  with  respect  to  the  underlying 
dominant  wind  waves. 

(iii)  surface  roughness  data  analysis 

As  a  component  of  our  data  analysis  effort,  we  continued  our  effort  to  develop  new  phase-resolved 
descriptions  of  sea  surface  roughness.  In  addition  to  state-of-the-art  Fourier  (e.g.  Elfouhaily  et  ah, 

2003)  and  wavelet  techniques,  we  are  investigating  various  ‘riding  wave’  analyses  (RWA).  This 
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approach  envisages  an  ‘individual  wave’  deeomposition  eapability  so  that  loeal  physieal  roughness 
elements  ean  be  deteeted  and  eharaeterized  along  with  their  spaee-time  phasing,  thereby  overeoming 
the  elassieal  Fourier  speetrum  issue  of  bound  versus  free  wave  eontributions  in  assessing  true  physieal 
sea  surfaee  roughness. 

IMPACT/APPLICATIONS 

This  effort  will  provide  a  far  more  detailed  eharaeterization  of  the  wind  driven  air-sea  interfaee, 
ineluding  wave  breaking  (whiteeaps  and  mieroseale  breaking).  This  is  needed  to  provide  more 
eomplete  parameterizations  of  these  proeesses,  whieh  will  improve  the  aeeuraey  of  oeean  optieal 
radiative  transfer  models  and  trans-interfaeial  image  reeonstruetion  teehniques. 

RELATED  PROJECTS 

The  present  projeet  is  related  generieally  to  our  eurrent  ONR  sea  surfaee  wave  projeet  in  the  CBLAST 
Hurrieane  DRI  entitled:  ‘Wave  breaking  influenee  in  a  eoupled  model  of  the  atmosphere-oeean  wave 
boundary  layers  under  very  high  wind  eonditions’.  While  the  wind  speed  regimes  in  RaDyO  and 
CBLAST  Hurrieanes  are  very  different,  eommon  elements  in  these  two  projeets  inelude  the  need  to 
better  understand  and  parameterize  the  breaking  proeess  and  how  it  oeeurs  at  the  different  wave  seales. 

Our  CBLAST  effort  has  resulted  in  a  eapability  for  foreeasting  wave  breaking  of  the  dominant  waves. 
These  foreeasts  validate  well  at  moderate  wind  speeds  of  around  12  m/s.  Validation  of  hurrieane 
breaking  waves  awaits  the  data  proeessing  by  other  CBLAST  Pis.  Our  effort  has  also  highlighted  the 
need  to  better  understand  breaking  at  the  shorter  seales,  where  the  breaker  frequency  statistics  appear 
to  fall  off  towards  shorter  seales.  However,  balaneing  wind  input  to  short  waves  with  breaker 
dissipation  rate,  present  modeling  suggests  that  short  wave  breaking  statisties  should  inerease  towards 
shorter  seales.  We  will  revisit  this  issue  with  the  new  insights  provided  by  our  RaDyO  datasets. 
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